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ABSTRACT 


r 

'YVETTE  is  a  free-fall  instrument  package  which 
measures  temperature,  conductivity,  pressure,  and  horizontal 
velocity  in  the  upper  ocean.  From  the  direct  measurements, 
the  density,  p,  and  gradient  (or  differenced)  quantities  such 
as  Brunt-Vaisala  frequency,  N=-(  g/p)  ( &p/.az) ,  and  the  square 
of  the  vertical  shear,  Sl2_(  az)^+(  6v/ az) can  be  cal¬ 
culated.  Data  from  YVETTE  have  been  used  to  determine 
empirical  relations  between  and  S^,  and  variations  in 
these  relations  among  different  oceanographic  regions. 

The  purpose  of  this  report  is  to  address  questions 
that  have  arisen  as  to  the  reliability  of  the  measurements 
and  the  ensuing  analyses.  Sources  of  noise  in  the  measure¬ 
ments  are  identified  and  the  noise  magnitudes  are  estimated. 
Estimates  of  noise  in  the  calculated  values  of  and  are 
made  based  on  an  expression  derived  by  Gregg  (1979). 
Relative  noise  errors  decrease  as  the  quantities  increase. 
As  increases  from  10"®  to  10"^  s"^,  the  relative  noise 
error  decreases  from  as  much  as  100%  to  as  little  as  2%. 
Similarly,  as  increases  over  the  same  range  of  values,  the 
relative  noise  error  decreases  from  30%  to  less  than  1%.  The 
apparently  large  noise  error  in  (compared  to  that  in  S^) 
is  due  to  the  low  precision  of  the  temperature  data  provided 
for  the  analyses.  -  It  is  concluded  that  the  shear  charac¬ 
terizations,  the  N^/S^  relations,  and  other  conclusions 
reached  to  date  are  not  affected  by  the  noise  errors  de¬ 
scribed  here.  Some  estimates  of  the  wavenumber  limits  of 
these  analyses  are  also  made. 
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Section  1 
INTRODUCTION 


1.1  BACKGROUND  AND  OBJECTIVES 

The  purpose  of  this  report  is  to  document 
identified  sources  of  error  in  data  from  the  free-falling 
density  and  velocity  profiler  YVETTE,  and  to  assess  the 
impact  of  these  errors  on  our  efforts  to  use  YVETTE  data  to 
characterize  upper  ocean  vertical  shear.  Vertical  profiles 
of  density  stratification  and  horizontal  velocity  obtained 
with  YVETTE  were  used  to  characterize  the  spatial  distribu¬ 
tion  of  small-scale  shear  in  the  upper  ocean,  to  investigate 
relationships  between  small  scale  vertical  shear  and  density 
stratification  (Patterson,  et  ,  1981a)  and  to  examine  the 
sensitivity  of  Richardson  nvmber  to  various  computational 
procedures  (Newman,  aJ* >  1981).  A  detailed  description  of 
YVETTE  is  presented  in  Evans,  £t  . >  (1979). 

Specific  objectives  of  this  study  were  to: 

•  estimate  the  minimum  vertical  scale  for  which 
YVETTE  conductivity  data  should  be  used  in  the 
computation  of  vertical  density  gradients, 

•  estimate  the  error  in  Br unt-Vaisala  frequency 
derived  from  YVETTE  measurements  of  temper¬ 
ature,  conductivity  and  pressure, 

•  estimate  the  specific  error  in  the  YVETTE 
acoustic  velocity  sensor,  and  the  overall  error 
in  computed  values  of  shear-squared  (S^), 
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•  assess  the  effects  on  calculated  of  the  loss 
of  a  fourth  decimal  place  in  the  original 
temperature  data.  (In  the  YVETTE  data  set 
supplied  to  SAI,  temperature  was  reported  to 
only  three  decimal  places.) 

1.2  MINIMUM  VERTICAL  SCALE  OF  YVETTE  CONDUCTIVITY 

DATA 


At  this  time  an  adequate  description  of  the 
actual  frequency  response  of  the  conductivity  cell  is  not 
available,  though  Michael  Gregg  (Applied  Physics  Laboratory/ 
University  of  Washington)  is  investigating  it.  Discussion 
of  this  objective  is  therefore  limited  to  the  following. 


A  statement  about  the  minimum  vertical  scale 
of  YVETTE  conductivity  data  requires  a  precise  knowledge 
of  the  frequency  response  of  the  sensor  and  the  fall  rate 
of  the  instrument.  For  an  ideal  instrument  with  a  flat 
frequency  response  in  the  pass  band  and  a  sharp  roll  off 

near  the  Nyquist  frequency,  the  data  are  useful  to  the 

vertical  scale  given  by  the  product  of  the  fall  rate  and 
the  reciprocal  of  the  Nyquist  frequency.  For  YVETTE  the 

sampling  frequency  is  2.5  Hz,  so  the  Nyquist  frequency 
is  1.25  Hz.  The  fall  rate  is  20  to  25  cm  s~  ^ .  If  the 
conductivity  response  were  ideal  as  described  above,  the 
minimum  vertical  scale  resolved  would  be  16  to  20  cm. 
Thus  1  m  resolution  would  be  within  the  capabilities  of 
YVETTE . 

De  facto  evidence  of  YVETTE's  ability  to  resolve 

small  scale  structure  in  the  conductivity  profile  exists. 
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Evans  (personal  communication)  notes  that  after  time-lag 
corrections  are  applied  to  the  temperature  profiles,  small 
scale  structure  (<  2  m)  observed  in  both  the  temperature  and 
conductivity  profiles  yield  calculated  density  profiles 
which  are  stable.  This  is  not  likely  to  be  fortuitous  and 
therefore  suggests  that  the  conductivity  sensor  on  YVETTE  is 
able  to  make  valid  measurements  at  vertical  scales  of  2  m. 

1.3  ERROR  IN  YVETTE  DATA 

The  overall  error  in  YVETTE  data  depends  on 
the  error  in  the  raw  data  introduced  by  the  imperfect 
response  of  the  sensor/recorder  package  as  well  as  the  way 
in  which  that  error  is  modified  by  subsequent  data  transfer 
and  processing.  Some  of  the  imperfect  response  can  be 
compensated  during  data  processing.  That  which  cannot  be 
compensated  contributes  to  measurement  error. 

Some  potential  sources  of  error  are  inherent 
within  the  sensor/recorder  package  before  it  is  deployed. 
These  include  zero  offset  or  drift,  calibration  error, 
sensor  inertia  (i.e.  time  response  and  sensitivity),  sensor 
precision,  and  least  bit  resolution  in  the  digital  recording 
system.  Other  sources  of  error  are  manifested  while  the 
instrument  is  collecting  data  and  are  the  result  of  extrane¬ 
ous  interactions  between  the  package  and  the  surrounding 
environment.  These  include  disturbance  of  the  environment 
by  YVETTE,  descent  attitude  problems  (such  as  a  mean  tilt  of 
the  body),  the  horizontal  accelerations  imposed  on  the  body 
of  YVETTE  by  the  large  scale  vertical  shear,  and  pendulum- 
type  body  motions.  Also,  the  descent  rate  decreases  some¬ 
what  as  the  density  of  the  surrounding  water  increases,  but 
this  has  a  negligible  effect  on  measurement  errors  discussed 
in  this  report. 
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For  the  computation  of  vertical  gradients  over 
relatively  small  vertical  scales,  the  important  consid¬ 
eration  is  not  absolute  accuracy  of  the  measured  quantities, 
but  the  ability  of  the  sensor  system  to  resolve  small  scale 
fluctvations  in  the  measured  quantities. 


Like  other  similar  profiling  devices  currently  in 
use  (cf.  Lambert  and  Patterson,  1980),  the  YVETTE  sensor 
package  is  still  under  development.  While  the  response 
characteristics  of  the  individual  sensors  are  reasonably  well 
established,  the  package  as  a  whole  has  undergone  only 
limited  testing.  There  are  some  aspects  of  the  response 
which  are  not  well  known  and  can,  at  best,  be  inferred  by 
indirect  methods. 


1.4  OVERVIEW 

Section  2  presents  an  estimate  of  the  measure¬ 
ment  errors  associated  with  the  conductivity,  temperature, 
and  pressure  sensors.  It  includes  a  discussion  of  the 
effects  of  various  data  processing  techniques,  and  an 
estimate  of  the  relative  error  in  derived  values. 

Section  3  presents  an  estimate  of  the  measure¬ 
ment  error  associated  with  the  acoustic  current  meter, 
a  discussion  of  the  effects  of  data  processing,  an  estimate 
of  relative  error  in  derived  S^  values,  and  a  discussion  of 
the  impact  of  these  errors  on  investigations  utilizing  YVETTE 
data. 


4. 


A  summary  and  conclusions  are  presented  in  Section 
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Section  2 


THE  EFFECTS  OF  MEASUREMENT  ERROR  IN  CONDUCTIVITY, 
TEMPERATURE,  AND  PRESSURE 

According  to  Gregg  (1979),  the  two  principal  error 
types  in  CTD/velocity  profilers  are  system  noise  and  dc 
offsets  (bias  errors).  System  noise  in  the  conductivity, 
temperature  and  pressure  measurements  introduces  relative 
noise  in  the  calculated  value  of  N^  of  0.4*  to  over  100" 
depending  upon  vertical  spacing  of  the  data  used  in  the 
calculation.  Dc  offsets  in  conductivity,  temperature  and 
pressure  introduce  relatively  small  errors  in  N^  of  0.5%  to 
4.5%  which  result  from  evaluating  the  nonlinear  equation  of 
state  at  incorrect  values  of  T,  C,  and  P. 

In  this  section,  Gregg's  derivation  for  the  effect 
of  system  noise  is  summarized.  His  equation  estimating  the 
relative  noise  in  N*  due  to  system  noise  is  presented  and 
applied  to  the  YVETTE  measurements.  It  is  shown  that  the 
relative  noise  in  N^  from  YVETTE  data  is  large  in  weak 
stratification  and  becomes  smaller  as  stratification 
increases . 

2.1  THE  CALCULATION  OF  ERROR  IN  N^ 

A  very  precise  method  for  calculating  the  Brunt 
Vaisala  frequency  (N)  from  data  is  given  (e.g.  Gregg,  1979) 
by  the  formula; 

ip  • 


i. 
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where  a  =  potential  specific  volume  referred  to  pressure 

midway  between  two  observations, 


Acc  *  change  in  potential  specific  volume  over  the 
sample  interval  using  as  reference  the  pres¬ 
sure  midway  between  the  two  observations, 


AP  =  difference  in  pressure  over  the  sample 
interval,  and 


g  =  acceleration  of  gravity. 

Gregg  uses  a  linearized  form  of  the  equation  of 
state  (relating  a  to  measured  values  of  temperature,  T,  con¬ 
ductivity,  C,  and  pressure,  P)  and  derives  an  expression  for 
,  the  root  mean  square  error  in  due  to  system  noise. 

The  procedure  involves  taking  the  ensemble  average  of 

individual  instances  of  error.  If  nis  is  the  error  intro- 

N 

duced  in  the  ith  instance  by  an  error  of  nl  in  Aa  and 

A  a 

Ti^p  in  AP,  then 


°N2 


1/2 


(2.2) 


From  equation 


^N2 


2.1: 

B  2 
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JP  + 

^Ip 

¥2 
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Gregg  first 


defines 
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(the  mean  square  error  in 
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2 

a 
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C  C  P  p 


(2.4) 
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where  2  =  mean  square  error  in  potential  specific 
volume , 


:^2  =  mean  square  error  in  observed  temperature 


=  mean  square  error  in  observed  conductivity 


:p2  s  mean  square  error  in  observed  pressure 


^  ^  J-  i:s 
rT  >S  3T 


3  j  S 

oS  ?c 


(2.5) 


—  _  -  r 

cS  :-P  ^  5T 


where 


S  =  salinity,  T=temperature,  C=conductivity  and 
P<=pressure . 


r  =  adiabatic  lapse  rate,  a  function  of  local 
T,  S  and  P. 


Then,  computing  and  from  and  np  and  substituting 

in  equations  2.4  and  2.3  gives; 


^2  2  +  -  2  2(AaJ  v  o  ^ 


a„2  Op2 


+  2  3 

(AP)2  (^)2 


2  41’ 

Op  (2.6) 


I 


where 


Aa  =  the  change  in  potential  specific  volume  as 

estimated  by  an  ensemble  average  of  observed 
changes,  and 

aP  =  the  change  in  pressure  as  estimated  by  an 
ensemble  average  of  observed  changes. 


The  last  two  terms  in  equation  2.6  are  second  order  and, 
according  to  Gregg,  can  be  neglected  without  changing 
estimates  of  c ^2  by  more  than  10%  if 


< 


0.1  or 


^  p  <  0.1. 

ZTp 


(2.7) 


Neglecting  these  second  order  terms,  dividing  by  n2  to 
provide  the  relative  error  of  and  substituting  the  system 
noise  contributions  for  the  specific  volume  noise,  Gregg 
derives  an  alternative  form  of  the  noise  equation, 


InH  -  ^ 

n2  “  AP 


(J) 


2  0^  + 
T  "T  ^ 


■(2.8). 


Though  the  coefficients  K-p,  Kc  and  Kp  are  each  functions 
of  T,  C  and  P,  the  coefficients  vary  only  slightly  over  the 
water  column  compared  to  variations  in  N.  Thus  Gregg  (1979) 
approximates  the  variable  values  of  the  coefficients  with 
constant  values,  though  he  does  not  indicate  his  choice. 
Gregg  does  provide  figures  that  illustrate  the  dependence  of 
Kx,  Kc  and  Kp  on  T,  S  and  P.  Assuming  typical  ocean  values 
of  T=10OC,  S=35°/^q  and  P=1000  dbar=107  Pa,  we  select 
the  following  values  for  the  coefficients  from  Gregg's 
figures: 
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r 


Kt  = 

8.7 

X 

10- 

m  ^ 

kg-l 

oc-  1 

Kc  = 

-7.0 

X 

10-6 

or 

m  ^ 

kg-i 

(S/m)-l 

-7.0 

X 

10- 

m  ^ 

kg-i 

(mmho/cm) 

Kp  = 

3.0 

X 

10-1“ 

or 

m  ^ 

kg-l 

Pa-  1 

3.0 

X 

10-10 

m  ^ 

kg-  1 

dbar-  1 

Gregg  provides  plots  of  the  relative  error  in 
as  a  function  of  N  and  AP  for  two  assumed  noise  levels 
in  CTD  data,  reproduced  here  as  Figures  2.1  (a)  and  (b) . 

2.2  System  Noise  in  YVETTE  CTD  Data 


The  raw  data  collected  by  YVETTE  consists  of 
triplets  of  temperature,  conductivity  and  pressure  measure¬ 
ments  recorded  at  a  rate  of  2.5  Hz.  Typical  fall  rates  are 
25  cm  s~  ^ ,  yielding  a  raw  data  triplet  every  10  cm.  The 
datum  from  each  measurement  is  recorded  in  a  16  bit  computer 
word,  i.e.,  the  digital  resolution  is  one  part  in  2^^  over 
the  range  of  the  sensors.  The  sensor  ranges  and  resulting 
parameter  resolutions  are  indicated  in  Table  2.1.  Also  given 
in  Table  2.1  are  values  of  quantizing  noise,  YVETTE  system 
noise,  and  noise  in  the  processed  data.  These  values  are 
derived  below. 

For  data  with  a  known  digital  resolution  the 
noise  can  be  no  less  than  a  value  called  the  quantizing 
noise,  oq.  It  is  given  by  the  relation: 

2 

2  (resolution ) 

Oq=  •  (2.9; 
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Table  2.1 

YVETTE  SENSCMl  SPECIFICATIONS 


Values  of  the  quantizing  noise  are  also  given  in  Table  2.1. 


Fofonoff  ^  (1974)  have  performed  spectral 
analyses  of  the  output  of  similar  temperature  and  conductiv¬ 
ity  sensors  from  a  Brown/WHOI  CTD  in  a  homogeneous  environ¬ 
ment  and  state  that  system  noise  Oq  associated  with  these 
sensors  is  equal  to  the  quantizing  noise  oq.  In  the  case 
of  temperature  from  a  prototype  CTD  they  did  find  system 
noise  greater,  by  about  a  factor  of  than  the  quantizing 
noise  but  they  further  state  that  the  problem  leading  to  the 
increased  noise  has  been  corrected.  Because  the  CTD  in 
YVETTE  is  a  very  early  version,  we  will  use  the  larger  value 
of  system  noise  (Vs*  x  oq)  for  temperature.  For  the  pres¬ 
sure  sensor  we  will  assume,  as  Gregg  (1979)  does,  that  the 
noise  is  0.029  dbar  which  is  higher  than  the  quantizing  noise 
level.  These  system  noise  (oq)  values  are  given  in  Table 
2.1. 


The  estimated  YVETTE  system  noise  levels  of  Table 
2.1  can  not  be  used  directly  in  either  equation  2.6  or  2.8 
because  the  raw  YVETTE  data  are  not  used  in  the  calcula¬ 
tion  of  N^.  The  temperature  data  are  first  corrected  for 
a  time  lag,  then  salinity  is  calculated  and  finally  the  T 
and  S  series  are  filtered  and  subsampled  at  1  m  (4  sec) 
intervals.  A  different  filter  is  applied  to  the  pressure 
series.  We  will  assume  that  calculating  salinity  has  no 
effect  on  the  manner  in  which  noise  in  temperature,  pressure 
and  conductivity  affects  noise  in  N^.  That  is,  we  will  per¬ 
form  this  error  analysis  as  if  T  were  time  lag  corrected, 
then  T  and  C  were  filtered,  then  P  was  filtered  and  then 
was  calculated  from  T,  C  and  P  series.  It  should  be  pointed 
out  that  filtering  the  data  reduces  the  signal  as  well  as  the 
noise  at  wavenumbers  above  the  filter  cutoff  wavenumber. 
This  is  an  unavoidable  result  of  filtering,  but  only  serves 
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to  set  a  high  wavenumber  limit  on  the  usable  data  and  is  not 
considered  a  source  of  error  in  this  discussion. 


The  time  lag  correction  increases  the  noise 
to  signal  ratio  in  a  manner  that  depends  upon  frequency 
(Fofonoff  et  ,  1974).  A  smoothing  (low  pass)  filter 

decreases  the  noise  to  signal  ratio  in  a  different  manner 
that  is  also  frequency  dependent.  The  frequency  responses 
of  these  operations  are  given  by  their  transfer  functions 
and  R2  respectively.  Assigning  system  noise  that  is 
originally  white  (i.e.,  frequency  independent)  with  original 
rms  noise  Oq  =  3x10“**  °C  (Table  2.1)  the  resulting  noise 
has  a  spectral  shape  Sr(w)  given  by 


Sj^((i))  =  Rj(a))  •  R2(  ti>)  . 


where  is  the  noise  per 

width  of  the  measurements, 
therefore  is  a  constant. 
( is  given  by 


unit  frequency 
The  noise  is 
The  resulting 


over  the  band 
assumed  white, 
noise  variance 


c 


R 


2 


R2(  u))  do). 


That  is,  the  change  in  variance  is  given  by  the  ratio 


f  R^(  (i))  •  R2(  w)  d  w. 
o 


(2.10) 


The  transfer  function  for  the  time  lag  correction,  Rj ,  also 
depends  upon  the  number  of  points  (m)  over  which  a  least 


squares  linear  fit  is  applied  in  order  to  estimate  the  time 
derivative  of  temperature.  According  to  Evans  (personal 
communication)  m=3  for  the  YVETTE  temperature  data.  liien, 
according  to  Fofonoff  ^  (1974): 

Rj  =  1  -  I  sin^  ^  ^  “l 


where 


C  TC  OJ 


wQ  =  Nyquist  frequency  of  time  series,  and 

nL  time  lag  in  number  of  time  intervals. 

Evans  f?'?*ls  values  of  hl  in  the  range  of  3.3  to  3.4  which 
for  a  sampling  interval  of  0.4  s  (i.e.,  a  sampling  fre¬ 
quency  itf.  V.5  Hz)  amounts  to  1.34  s.  It  must  be  pointed 
out  that  this  is  several  times  larger  than  the  lags  of 
~0.16  s  determined  for  similar  instruments  (e.g.  by  Fofonoff 
et  al . ,  1974,  and  Joyce,  1976).  Such  lags  were  deemed  ne¬ 
cessary  and  reasonable  in  order  to  eliminate  density  inver¬ 
sions  in  the  calculated  density  profile  (Evans,  personal 
communication) . 

Evans  smoothed  the  YVETTE  data  with  a  time  domain 
filter  having  a  frequency  response  similar  to  the  four  pole 
Butterworth  filter  whose  transfer  function  is 


1  +  (20  6  A)’ 


We  numerically  integrated  the  product  of  Rj  and 
R2  over  the  interval  0_<a)^(»p(0£6£  11/2)  and  obtained 
for  the  ratio 


2-10 


0. 3865,  or  a 


0.622 


Thus,  using  Oq  for  temperature  from  Table  2.1,  the 
resultant  rms  noise  in  temperature  is  1.9xl0~‘*oc  in  the 
smoothed  1  m  temperature  series. 


The  conductivity  data  were  not  lag  corrected 
so  Ri  =  1  for  all  w  and 

2 

-Sjr  =  0.1612, or  Oq  =  0.401  a  . 

£t  tl  O 

a 

o 


Thus,  referring  to  Table  2.1  the  resultant  rms  noise  in  con¬ 
ductivity  is  1.2x10“**  mmhos  cm"  ^  in  the  smoothed  1  m 
conductivity  series. 


The  pressure  was  smoothed  by  applying  a  99  point 
filter  to  the  time  series.  The  weights  of  the  filter  were 
obtained  from  the  Gaussian  distribution  for  the  interval  -2d 
to  +2d ,  where  d  represents  the  standard  deviation  of  the 
distribution.  For  99  points  to  equal  4d  implies  that  d 
equals  9.8  seconds,  or  24.5  sampling  intervals.  Taking  the 
Fourier  transform  of  the  Gaussian  and  squaring  it  gives  the 
response  function  (R3)  of  this  filter: 


R 


2  2 
b) 


As  before  we  integrate  this  over  the  interval  to 

obtain  for  pressure: 
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4 


R 

— =•  =  .0058  or  Oq  =  .076  o 

21  ri  O 

Referring  to  Table  2.1,  Oq  for  pressure  is 
2.9x10“^  dbar  so  according  to  the  above  result,  the  resultant 
noise  in  the  processed  data  is  2.2x10“^  dbar. 

The  system  noise  estimates  for  the  processed  data 
are  given  in  the  final  column  of  Table  2.1.  These  noise 

levels  can  then  be  substituted  in  equation  2.8.  We  can  apply 
equation  2.8  if  the  condition  of  equation  2.7  is  met.  YVETTE 
meets  the  pressure  condition  if  aP  >_  10  V  2  op,  i.e, 
AP  2  0.41  dbar.  So,  for  separations  in  excess  of  0.41  dbar 
equation  2.8  can  be  used.  'Hie  1  m  data  obviously  satisfy 
this  requirement.  The  relative  error  in  as  a  function  of 
N  and  AP  is  plotted  for  the  processed  YVETTE  data  in  Figure 
2.2.  For  comparison  a  curve  from  Figure  2.1(a)  is  also 
plotted  showing  that  the  processed  YVETTE  data  does  not  per¬ 
form  as  well  as  unsmoothed  CTD  data  with  only  quantizing 
noise.  One  dbar  separations  give  errors  less  than  10%  in 
onlf  in  strongly  stratified  regions  such  as  the  seasonal 
thermocline.  In  the  less  stratified  main  thermocline  rela¬ 
tive  errors  of  less  than  20%  occur  only  for  separations  of  4 
dbar  or  greater  in  N^.  In  relatively  homogeneous  regions 
like  the  deep  ocean  or  mixed  layers,  even  at  4  dbar  separa¬ 
tions,  relative  errors  can  be  as  large  as  70%. 

2.3  ERROR  IN  YVETTE  DATA  PROVIDED  TO  SAI 

The  data  provided  to  SAI  for  analysis  are  not  iden¬ 
tical  to  Evans'  processed  YVETTE  data.  Temperatures  were  not 
provided  to  10"‘*oc  to  take  full  advantage  of  the  low  noise 
level  (~1 .9x10“  **00  indicated  in  Table  2.1.  Instead,  the 
temperatures  were  provided  to  the  nearest  10“ ^oc  which 
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Figure  2.2  Relative  error  in  as  a  function  of  N. 

Solid  lines:  YVETTE  data  after  lag 
corrections  and  smoothing.  Three 
different  differencing  intervals  are  shown 


Dashed  line:  unsmoothed  CTD  data  with 
only  quantizing  noise  contributing. 


increases  the  quantizing  noise  to  2.9xl0“‘*oc.  This  inde¬ 
pendent  source  increases  the  total  rms  noise  of  the  data. 

6T0TAL  =  =  3.5  x  lO'^^C. 

In  addition,  salinity  resolved  to  1  ppm  was  provided  instead 
of  raw  conductivi ty  resolved  to  .001  mmhos  cm"^.  This  latter 
substitution  should  not  increase  the  noise  in  the  calculated 
values  of  a.  Pressure  values  were  reported  to  the  nearest 
0.1  dbar  which  implies  a  quantizing  noise  of  .029  dbar.  This 
substantially  degrades  the  pressure  series  from  the  level 
that  appears  in  Table  2.1. 

Substituting  the  larger  rms  temperature  and 
pressure  noise  values  in  equation  2.8  yields  a  relative 
error  of  (as  a  function  of  N  and  aP)  approximately  1.6 
times  higher  than  in  Figure  2.2  for  N<.01  rad  s~ ^  (5.5  cph) 
and  about  10  times  higher  for  N>.01  rad  s~  ^ .  This  case 
is  plotted  in  Figure  2.3.  For  values  of  N  less  than  5x10“^ 
rad  s~  ^  (3  cph),  a  4  dbar  spacing  is  required  to  ensure 
relative  errors  of  less  than  10%  in  .  In  regions  of  lower 
stratification  much  larger  separations  are  required. 

In  conclusion,  the  YVETTE  data  collection  and 
processing  procedures  do  not  compare  favorably  to  optimum 
CTD  performance  except  when  N  is  large.  Using  more  modern 
temperature  sensor  electronics  will  greatly  improve  the 
calculation  of  .  It  is  undesirable  to  round  temperature 
values  to  the  nearest  millidegree  or  pressure  values  to  the 
nearest  tenth  decibar,  since  this  step  further  degrades  the 
calculation  of  N  at  all  values  of  N. 
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N  (radians /second) 


Figure  2.3  Same  as  Fig.  2.2  except  that  the  solid  lines 
represent  the  performance  of  data  provided  to 
SAI  on  digital  tape.  The  increased  noise  level 
in  n2  is  due  to  increased  noise  in  the  tempera¬ 
ture  data  due  to  quantizing  noise  introduced  when 
temperature  values  are  rounded  to  the  nearest 
millidegree,  and  pressure  to  the  nearest  tenth 
of  a  decibar. 


Section  3 

THE  EFFECTS  OF  MEASUREMENT  ERROR  IN  HORIZONTAL  VELOCITY 


3.1  ERROR  ASSOCIATED  WITH  THE  VELOCITY  SENSOR 


The  YVETTE  acoustic  current  meter  (Figure  3.1)  was 
developed  by  T.  Gytre  at  the  Christian  Michelsen  Institute 
in  Bergen,  Norway.  It  consists  of  two  orthogonal  pairs  of 
ultrasonic  (~4MHz)  pingers.  Each  pair  measures  one  compo¬ 
nent  of  horizontal  water  velocity  relative  to  the  body  of 
YVETTE.  For  a  given  pair  of  transducers,  pulses  emitted 
simultaneously  from  each  pinger  propagate  in  opposite 
directions  along  the  straight  line  path  between  them.  The 
difference  in  travel  time,  t  ,  is  related  to  the  component 
of  mean  velocity  parallel  to  the  propagation  path  by 


_  2LV  cos  6 

r2  ^2  2  , 

c  -V  cos  r 


(3.1) 


where  L  is  the  separation  distance  between  pingers  (  30  cm), 
V  cos  f  is  the  component  of  mean  water  velocity  parallel  to 
the  propagation  path,  and  C  is  the  loc<il  speed  of  sound. 
Since  V  cos-  is  much  less  than  C,  eouution  (3.{)  can  be 
approximated  by 


_  —  2LV  cos  6 


(Evans  et  al.,  1979). 


(3.2) 
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In  order  to  resolve  current  speeds  as  small  as  1  mm  s~l , 
the  noise  level  of  the  current  meter  should  be  no  more  than 
0.5  mm  s“l .  This  requires  a  resolution  of  Z,  t  of  10"10  s. 
In  order  to  achieve  this  high  degree  of  resolution,  the 
circuitry  has  been  designed  to  continually  monitor  and 
compensate  for  drift  in  the  electronic  zero  and  for  the 
variation  with  depth  of  the  speed  of  sound. 

The  travel  time  differences  are  measured  by 
letting  the  received  signal  from  one  pinger  start  a  high 
speed  ramp  and  the  signal  from  the  opposite  pinger  stop  it. 
The  ramp  amplitude  is  modified  by  the  local  value  of  sound 
speed  which  is  determined  from  the  total  travel  time  of  the 
pulse  between  pingers.  This  time  interval  is  measured  by  an 
oscillator-controlled  counter  which  is  started  when  the 
pulse  is  generated  and  stopped  when  the  pulse  is  received  at 
the  opposite  pinger.  The  least  bit  resolution  for  total 
travel  time  corresponds  to  a  least  bit  resolution  in 
sound  speed  of  approximately  11  m  s“l(-.75%)  or  a  quantiza¬ 
tion  noise  level  (i.e.  (least  bit)  of  approximately 
3  m  s~l(~.2Si).  From  equation  (3.2)  this  noise  level  in  the 
sound  speed  measurement  will  contribute  a  noise  of~.4%  to 
the  measured  water  velocities.  Measured  water  velocities 
are  typically  of  order  1  cm  s~l.  This  is  because  YVETTE 
is  advected  horizontally  by  the  large  scale  velocity 
structure  and  that  which  is  actually  measured  is  the  water 
velocity  relative  to  the  body  of  the  instrument.  In  a  worst 
case  situation  where  a  large  relative  velocity  of  say 
10  cm  s~l  is  measured,  the  error  due  to  the  quantization 
noise  in  the  sound  speed  determination  is  only  0.04  cm  s“l. 
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■nils  and  any  other  source  of  noise  in  the  in¬ 
dividual  velocity  measurements  are  reduced  by  the  averaging 
that  occurs  when  the  data  is  recorded.  While  velocities  are 
sampled  at  a  rate  of  30  Hz,  they  are  digitally  recorded 
with  a  least  bit  resolution  of  0.05  cm  s~  ^  at  a  rate  of 
2.5  Hz.  Between  recordings,  travel  time  differences  are 
accumulated  and  that  which  is  actually  recorded  is  an  average 
velocity  over  12  observations.  Since  YVETTE  descends  at  a 
rate  of  approximately  25  cm  s~  ^ ,  this  corresponds  to  an 
average  over  a  vertical  distance  of  approximately  10  cm.  The 
noise  level  associated  with  the  average  value  is  therefore 
smaller  than  that  for  any  individual  value  by  a  factor  of 
l/VlS  and  becomes  1.4x10"^  cm  s~^. 

The  net  noise  level  in  the  relative  velocity 
data  recorded  by  YVETTE  is  approximately  0.05  cm  s~  ^ 
(Evans  ^  ji*  »  1979).  This  is  equivalent  to  the  least  bit 
resolution  of  the  digitally  recorded  data  (Table  4.1).  The 
absolute  accuracy  is  not  as  good  as  this. 

Calibration  checks  on  the  accuracy  of  YVETTE 
velocity  measurement  were  attempted  in  a  tow  tank  at  the 
University  of  Rhode  Island.  The  plan  was  to  suspend 
YVETTE  from  a  trolley  and  tow  the  instrument  at  a  known 
constant  velocity  over  a  distance  of  15  m.  Unfortunately, 
the  tow  speed  was  variable  by  +  10%  thereby  necessitating 
a  very  large  number  of  runs  so  that  a  reliable  mean  value 
could  be  obtained.  The  proportionality  constant  between 
actual  and  measured  speed  differed  by  9%  from  the  value 
which  had  been  established  during  calibration  in  Norway 
two  years  earlier.  However,  because  of  potential  inade¬ 
quacies  of  the  calibration  facility  at  Rhode  Island, 
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David  Evans  (personal  communication)  elected  to  use  the 
Norwegian  calibration  results.  Hence,  the  accuracy  of 
the  measured  relative  velocities  is  not  known  exactly,  but 
is  believed  to  be  better  than  10%. 


A  calibration  error  is  a  systematic  error  and,  as 
such,  it  will  modify  in  a  linear  fashion  all  observed 
velocities  and  thereby  will  modify  all  values  of  and  Ri  . 
It  will  not,  however,  affect  the  qualitative  behavior  of 
these  parameters.  In  particular,  it  will  not  affect  rankings 
of  various  oceanic  regimes  or  areas  with  respect  to  levels  of 
^  or  Ri  (i.e.  the  calibration  error  will  cause  all  levels  to 
be  either  too  high  or  too  low).  It  will  also  not  affect  the 
correlation  between  N^  and  S^. 

The  orientation  of  the  body  of  YVETTE  with  respect 
to  the  earth's  magnetic  field  is  sensed  with  a  compass  which 
consists  of  a  bar  magnet  pivoted  on  fixed  bearings.  The 
angle  is  measured  capacitively  and  is  stored  digitally  with 
a  least  bit  resolution  of  approximately  3°.  Since  YVETTE 
rotates  very  slowly  with  depth  (approximately  1  rotation  per 
300  m) ,  compass  inertia  is  not  believed  to  be  a  problem. 
The  overall  noise  level  in  direction  is  not  known,  but  it  is 
probably  dominated  by  least  bit  resolution  (Table  4.1). 

3.2  ERROR  INDUCED  BY  INTERACTIONS  BETWEEN  YVETTE 

AND  THE  ENVIRONMENT 

YVETTE  has  been  designed  to  minimize  disturbance 
of  the  environment  at  the  point  where  observations  are 
obtained.  The  acoustic  transducers  are  approximately  60  cm 


3-5 


below  the  pressure  case  which  constitutes  the  main  body  of 
YVETTE  (Figure  3.1).  The  transducer  pairs  are  separated 
horizontally  by  approximately  30  cm  or  about  twice  the 
diameter  of  the  pressure  case.  Since  the  sensors  are  at  the 
lower  end  of  YVETTE,  they  penetrate  undisturbed  water  well 
ahead  of  the  main  body. 

Because  YVETTE  descends  at  a  relatively  slow  rate 
(~25  cm  s~^)  it  experiences  horizontal  accelerations  imposed 
by  the  local  velocity  field.  YVETTE  responds  as  a  high-pass 
filter  and  is  unable  to  directly  resolve  fluctuations  in 
horizontal  velocity  over  vertical  scales  greater  than  20  m. 
This  cutoff  is  established  by  examining  vertical  w<2venumber 
spectra  of  velocity  profiles  measured  with  YVETTE  (Figure 
3.2).  For  wavenumbers  greater  than  0.05  cycles  per  meter 
(i.e.  wavelengths  less  than  20  m)  the  spectrum  exhibits  the 
expected  slope  of  -2  to  -2.5.  However,  for  wavenumbers  less 
than  0.05  cycles  per  meter,  where  the  instrument  should  act 
as  a  high-pass  filter,  the  spectrum  is  indeed  flat.  Thus  for 
wavenumbers  higher  than  0.05  m“  ^ ,  YVETTE  directly  measures 
velocity  subject  to  the  rather  small  errors  discussed  in  this 
Section.  For  wavenumbers  lower  than  0.05  m" ^ ,  an  algorithm 
is  required  to  reconstruct  the  long  wavelength  components  of 
the  velocity  profile.  David  Evans  has  constructed  such  an 
algorithm  which  reconstructs  the  actual  velocity  profile  in 
an  incremental  fashion  by  predicting  the  response  of  YVETTE 
to  the  sequence  of  observed  relative  velocities.  This 
algorithm  is  described  in  Rubenstein,  »  (1981).  It  has 
proven  to  be  capable  of  reproducing  to  within  2-3  cm  s" ^  rms 
the  velocity  structure  recorded  by  a  simultaneous  Absolute 
Velocity  Profiler  (AVP)  drop  performed  by  T.  Sanford.  The 
rms  velocity  difference  between  AVP  and  corrected  YVETTE 
profiles  is  due  to  the  high  wavenumber  content  of  the  YVETTE 


3-6 


profile  (which  is  lacking  in  the  lower  resolution  AVP 
profile)  ,  to  the  error  at  low  wavenumbers  in  the  corrected 
YVETTE  profile,  and,  of  course,  to  measurement  noise  in  the 
AVP  and  YVETTE  profiles.  According  to  Evans  (private 
communication,  1982),  much  of  the  rms  velocity  difference 
between  the  AVP  and  corrected  YVETTE  profiles  is  due  to  the 
high  wavenumber  content  of  the  YVETTE  profile.  Because  the 
error  in  the  corrected  YVETTE  velocity  profile  is  not  known 
for  the  band  of  wavenumbers  less  than  0.05  m~^,  but  is 
thought  to  be  small,  we  assume  this  error  to  be  negligible. 

According  to  Evans  ^  al.  (1979),  the  expected  mean 
value  of  either  velocity  component  observed  by  YVETTE  is 
zero.  This  is  due  to  the  fact  that  YVETTE  tends  to  rotate 
slowly  as  it  descends  and  also  tends  to  high-pass  filter  the 
actual  velocity  profile  to  scales  smaller  than  the  rotation 
scale.  Calculations  of  this  mean  over  many  segments  of  many 
profiles  have  shown  that  it  tends  to  vary  about  zero  by 
about  0.1  cm  s“^.  The  reason  for  this  apparent  drift  is  not 
clear.  According  to  Evans  (personal  communication)  this 
offset  varies  very  slowly  with  depth.  It  therefore  has  a 
long  vertical  wavelength  and  cannot  contribute  significantly 
to  the  noise  level  in  shear  over  small  scales. 

The  possibility  of  momentary  instrument  tilt  due 
to  strong  local  shear  was  considered  by  Evans,  ^  (1979). 
They  calculate  a  tilt  of  0.03°  in  response  to  a  rather 
strong  shear  of  5  x  lO”^  s"^.  A  tilt  of  this  magnitude  would 
induce  a  velocity  error  of  1.25  x  10"^  s“ ^  in  the  direction 
of  the  shear.  A  shear  of  5  x  10“^  s~  ^  corresponds  to  a 
velocity  difference  of  10  cm  s"^  over  half  the  Instrument 
length  (2  meters).  The  relative  velocity  measured  by  YVETTE 
would  therefore  be  10  cm  s“ ^  so  the  error  due  to  tilt  in  this 
extreme  case  is  of  order  0.1%. 


If,  while  descending,  YVETTE  were  to  have  a  mean 
tilt  due  to  asymmetries  in  the  distribution  of  mass  within 
the  instrument  package,  this  would  cause  an  apparent  constant 
zero  offset  in  one  or  both  velocity  components.  At  a  descent 
rate  of  25  cm  s“  ^  the  0.1  cm  s“  ^  value  cited  above  would 
correspond  to  a  mean  tilt  of  about  0.2°.  To  eliminate  such 
effects,  the  mean  for  each  profile  is  removed  from  the  data 
before  other  processing. 

From  vertical  wavenumber  spectra  of  YVETTE  velocity 
profiles  (Figure  3.2),  it  has  been  deduced  that,  while  de¬ 
scending,  the  tube  sometimes  exhibits  pendulum-like  oscil¬ 
latory  motions  with  a  period  of  approximately  6  seconds. 
This  type  of  response  is  apparently  common  for  long  cylindri¬ 
cal  tubes.  In  a  worst  case  these  oscillations  were  mani¬ 
fested  in  the  velocity  profile  with  an  amplitude  of  0.25  cm 
s~^.  A  digital  notch  filter  has  been  developed  and  used 
successfully  in  removing  this  effect  from  the  data.  In  the 
data  set  used  for  our  analyses  this  oscillation  was  removed 
by  smoothing  the  data  over  a  vertical  scale  of  2  m.  A  de¬ 
scription  of  this  smoothing  procedure  is  presented  in  Lambert 
et  ( 1980) . 

The  identifiable  sources  of  noise  due  to  inter¬ 
action  of  YVETTE  with  the  environment,  have  been  made  insig¬ 
nificant  either  through  instrument  design  or  through  sub¬ 
sequent  data  processing.  Shear  spectra  show  the  expected 
rise  with  +2  slope  at  high  wavenumbers  only  when  the  signal 
energy  level  drops  below  that  associated  with  least  bit 
noise  of  +  0.05  cm  s“  ^  in  the  velocity  measurements 
(Evans  e^  1979).  It  can  therefore  be  concluded 
that  YVETTE  can  resolve  vertical  shear  over  vertical 
scales  of  2  to  20  m  with  a  noise  level  which  is  consistent 
with  the  least  bit  resolution  in  velocity.  The  noise  level 
in  shear  and  are  derived  in  the  next  section.  Tbe 
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absolute  accuracy  of  shear  over  these  small  scales  may  have 
a  systematic  error  as  large  as  10%  due  to  a  possible  error 
in  the  calibration  constant.  Over  larger  vertical  scales 
both  precision  and  accuracy  of  the  shear  values  are  reduced 
because  of  the  horizontal  drag  on  YVETTE  which  can  be  only 
approximately  modeled  and  because  of  the  apparent  zero 
drift  which  is  due  to  unknown  causes. 

3.3  ESTIMATES  OF  ERROR  IN  S2 

The  error  in  S2  over  small  scales  derives  from 
the  noise  level  of  the  velocity  sensor  and  a  possible 
systematic  calibration  error  as  discussed  above.  In  the 
data  set  supplied  to  SAI,  velocity  was  reported  to  10“2  cm 
s~l.  Since  this  precision  is  better  than  the  0.05  cm  s-1 
least  bit  resolution,  it  does  not  contribute  to  error  in 
velocity  or  $2. 

The  Butterworth  smoothing  filter  that  was  applied 
to  the  profiles  effectively  removes  fluctuations  with 
vertical  scales  less  than  2  m  and  reduces  the  noise  level  in 
the  profile  by  a  factor  of  0.401  (cf.  discussion  on  the 
effects  of  the  Butterworth  filter  in  Section  2.2).  Hence 
the  rms  noise  level  in  the  filtered  velocity  data  is 

Oy  =  <Eu>^  =  0.401  X  0.05  =  2.0  x  10*2  cm  s-1 

where  Ey  represents  the  noise  error  in  one  sample  and 
the  brackets  denote  an  ensemble  average  over  many  samples. 
The  components  of  vertical  shear  are  computed  by  simple 
differences,  hence 

iii  ■  ^  -  u(z  +  Az)  -  u(z) 

Sz  Az  ~  Az 
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The  noise  error  associated  with  -u  is  Oy  V?,  assumiftr^  the 
noise  in  u(z  +  Az)  to  be  uncorrelated  with  the  noise  irs  u(z). 
The  noise  level  in  one  component  of  vertical  shear  is  thus 


V3 0  2.8x10' 

IS  ""  "  - ir¬ 


es. 3) 


where  A  z  is  expressed  in  meters. 

For  A  z  =  2  m, 

=  1.4  X  10-4  s-1  . 

Vertical  shear  squared,  S2,  is  given  by 
S2  =  u|  +  v|  . 

For  purposes  of  this  error  analysis  it  is  assumed  that 

<U2>"  =  <v|>^  or  S2  s  2u^  ,  (3.4) 

and  in  the  remaining  discussion  both  components  of  shear 
will  be  represented  by  Uz> 


Now  according  to  Hajek  (1969)  a  general  function 
g(x),  of  a  random  variable  has  the  following  expectation 
value u  and  variance  c 2; 

t(g(x))  ss  g(u(x))  (3.5a) 

2 

a^(g(x))  as  [g’(V'(x))]  a^(x)  (3.5b) 

provided  that 


g'(x)  «  constant  for  jx  -  u(x)  |  <  3c(x)  .  (3.5c) 
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Now  let 


g(x)  =  S2  (EuJ 


where  according  to  (3.4) 


S2(Eu  )  =  2(u2  +  E„  )‘ 


Then  g'(x)  becomes 


(3.6) 


_d 

dE. 


u. 


("uj 


so  (3.5c)  is  satisfied  provided  Eu^<<Uz.  From  (3.3)  and  for 
Lz  s  2m  this  condition  becomes  U2;>>1.4  x  10“'^  s“l.  Inspec¬ 
tion  of  the  YVETTE  data  shows  this  condition  to  be  satis¬ 
fied  in  nearly  all  observations.  The  rms  noise  in  $2  can 
therefore  be  estimated  using  (3.5b).  Substituting  (3.6) 
into  (3.5b)  the  mean  square  noise  is  then 


(3.7) 


But 


u(Eu  )  =  0 
z 


and 


.  (^\ 


I  “  ■ 


So 


•  7^]  'u  - 


(3.8) 
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Substituting  (3.4)  into  (3.8) 


43 

iz 


a 


u 


or 


8  X  10 


-4  ... 


_z 


(3.9) 


Hence  the  noise  level  in  S2  is  a  function  of  both  the 
vertical  differencing  interval,  -z,  and  the  magnitude  of  the 
shear  itself.  For  Lz  =  2m, 


^g2  *  (4  X  10“4  s“^)S  . 

The  way  in  which  the  relative  noise  level,  ^s2/s2,  varies 
with  $2  is  illustrated  in  Figure  3.4.  For  Lz  =  2  m, 
the  noise  level  exceeds  the  signal  level  for  values  of 
$2  less  than  1.0  x  lO"^  s”2  (or  values  of  S  less  than 
4.n  X  10-4  s-1). 


The  accuracy  of  each  component  of  vertical  shear 
over  small  scales  may  be  in  error  by  as  much  as  9%  due  to  a 
possible  calibration  error  in  the  YVETTE  velocity  sensor. 
However,  as  noted  earlier,  this  would  be  a  systematic  error 
which  would  cause  all  shear  values  to  be  consistently  high 
or  low.  A  9%  error  in  each  component  of  vertical  shear 
would  lead  to  approximately  18%  error  in  S2. 
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Section  4 

SUMMARY  AND  CONCLUSIONS 


This  report  has  addressed  questions  regarding  the 
magnitude  of  noise  in  YVETTE  measurements  of  temperature, 
conductivity,  pressure,  and  velocity  and  resulting  noise  in 
calculations  of  the  square  of  Br unt-Vaisala  frequency  (N^) 
and  shear  magnitude  squared  (S^).  Table  4.1  summarizes  the 
noise  estimates  for  each  of  the  measured  parameters.  Impli¬ 
cations  of  the  estimated  noise  in  and  for  recent 
analyses  are  discussed  later  in  this  section. 


The  combined  noise  from  YVETTE  temperature,  con- 


ductivitv,  and  pressure  signals  produces  a  relative  noise 


sample  spacing  corresponding  to  a  pressure  difference  of  AP  = 


2  dbar: 


(N<2xl0‘^rad 


-2^  =  1.57x10"^,  (N>4xl0"^rad  s~^  ]  ■ 


For  values  of  N  ^etween  2x10“^  rad  s~  ^  and  4x10“**  rad  s~  ^  the 
dependence  of  ^  undergoes  a  smooth  transition  from  varying 
as  to  being  independent  of  N.  For  vertical  sample 
spacings  corresponding  to  pressure  differences  other  than  2 
dbar,  the  above  relative  noise  values  must  be  multiplied  by 


(■  2dbar  ] 
AP 
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In  the  YVETTE  data  provided  to  SAI  for  analysis, 
temperature  was  not  reported  to  sufficient  precision.  This 
effectively  increased  the  quantization  noise  for  temperature 
and  increased  the  relative  noise  as  follows  for  AP  =  2  dbar: 


2.2xl0~^frad  s~^  1 


fN<10~^rad  s"^  ), 


^  =  2x10 


(N>2xl0~^rad  s"^ ). 


As  in  the  previous  case,  the  dependence  of  relative  noise  ^ 
on  N  undergoes  a  smooth  transition  over  the  range  of  N  be¬ 
tween  10~2  and  2x10"^  rad  s“^.  For  other  vertical  spacings 

the  relative  noise  values  must  be  adjusted  by  ( 2dbar ) 

LP 

as  before. 


An  evaluation  of  several  types  of  velocity  measure¬ 
ment  error  concluded  that  the  only  important  sources  of  error 
are  calibration  errors  of  unknown  magnitude  and  instrumental 
noise  which  is  dominated  by  least  bit  quantization  in  the 
digital  recording  system.  It  is  worth  noting  that  the  error 
due  to  application  of  Evans'  motion  correction  algorithm  has 
not  been  determined,  but  is  assumed  to  be  quite  small.  This 
assimption  is  based  on  a  comparison  of  a  profile  obtained 
from  Sanford's  Absolute  Velocity  Profiler  and  a  simultaneous 
YVETTE  profile  to  which  Evans'  motion  correction  algorithm 
was  applied. 

The  total  error  in  over  small  scales  derives 
from  the  noise  level  of  the  velocity  sensor  and  a  possible 
calibration  error.  The  noise  level  in  is  estimated  to  be 
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]  s  . 


2  rSxlO'^m  s 

°s  ‘  ^ - li — 

The  relative  noise  level  therefore  becomes  greater  than 

o 

1.0  for  values  of  less  than  1.6xl0~"  s“ ^  (or  values  of  S 
less  than  4.0x10“'*  s“^). 

Results  of  recent • efforts  to  calibrate  YVETTE  velo~ 
city  sensors  differed  by  approximately  9%  from  calibration 
results  two  years  earlier.  Because  the  more  recent  calibra¬ 
tion  was  less  reliable,  the  earlier  calibration  results  have 
been  used  in  processing  all  YVETTE  data.  A  9%  error  could 
cause  an  18%  shift  in  the  mean  and  variance  of  histograms  of 
and  a  shift  in  statistics  of  Richardson  number,  but  it 
would  not  change  the  form  of  observed  distribution 
functions. 


While  measurement  errors  can  affect  quantitative 
conclusions  resulting  from  analyses  of  YVETTE  data,  the 
qualitative  nature  of  the  conclusions  of  the  SAI  character¬ 
izations  suggest  that  they  are  not  significantly  affected. 

In  Patterson,  ^  (1981a  and  1981b),  the  rms 

error  in  for  a  vertical  difference  interval  Az  *  2  m  was 
estimated  to  be  1.75x10”^  s“ ^  which  is  slightly  lower  than 
2.38x10“®  s“  ^  from  Figure  2.3  for  N  <  10“^  s“  ^ .  For  N  >  10“^ 
s“^  (n2  >  10“**  s“^)  the  error  is  estimated  at  a  fixed  2%  of 
.  The  rms  error  estimated  for  ( Az  ■  2m)  in  Patterson, 
et  al .  (1981a  and  1981b)  was  3.4x10“®  s” ^  while  this  report 
estimates  the  error  to  be  (4x10“**  s“^)S  increasing  linearly 
with  S.  These  differences  in  estimated  error  are  not  large 
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enough  to  significantly  affect  the  results  and  conclusions  of 
Patterson,  ^  (1981a  and  1981b). 

In  Newman,  et  (1981),  the  rms  error  in  for 
Az  =  2  m  was  estimated  to  be  1.35x10" ®  s~ ^  and  the  rms  error 
in  to  be  (2.24x10"^*  s"  ^)S  +  1.25x10"®  s"^.  Again,  these 
estimates  are  approximately  a  factor  of  two  smaller  than  our 
current  estimates,  but  the  results  and  conclusions  of  that 
report  are  essentially  unchanged. 
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